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jrt ' Abstract. A veritable zoo of different knots is seen in the ensemble of looped 

-ij , polymer chains, whether created computationally or observed in vitro. At short 

^^ , loop lengths, the spectrum of knots is dominated by the trivial knot (unknot). The 

fractional abundance of this topological state in the ensemble of all conformations 

of the loop of A'' segments follows a decaying exponential form, ~ exp {—N/Nq), 

where Nq marks the crossover from a mostly unknotted (ie topologically simple) 

to a mostly knotted (ie topologically complex) ensemble. In the present work we 

'"^ , use computational simulation to look closer into the variation of A^o for a variety 

&^ ■ of polymer models. Among models examined, Nq is smallest (about 240) for 

Q ' the model with all segments of the same length, it is somewhat larger (305) for 

fi ' Gaussian distributed segments, and can be very large (up to many thousands) 

when the segment length distribution has a fat power law tail. 
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1. Introduction: Formulation of the Problem 

Of interest to anglers seeking to fill their creels and children seeking to fasten their 
shoes, a wide audience has found knots compelling from time immemorial. In 
the scientific community, knots have been featured in initial formulations of the 
nature of atoms, P (see a popular historical account in [2]), the formulation of 
certain path integrals, [3], and also in quantitative biology, where knots have been 
observed in, @IIH], and tied into, DNA, [HIIZI, where the space of knots is biologically 
created and manipulated [H|- Knots also have been observed occasionally in proteins, 

ma nn [12 ESI El]. 

Historically, the classification of knots and study of knot invariants were the 
first subjects of knot theory [2], and this remains in the center of attention among 
knot theorists of mathematical orientation Pj. Another fundamental aspect of knot 
theory is that of knot entropy. Physically, this group of problems comes to the 
fore in the context of polymers and biophysics. Mathematically, this issue belongs 
to both topology and probability theory and seems to remain underappreciated in 
the mathematics and mathematical physics community. Even the simplest question 
in this area is poorly understood: what is the probability that a randomly closed 
loop in 31? will be topologically equivalent to plane circle? In other words, using 
professional parlance of the field, what is the probability that random loop is a trivial 
knot (unknot), Oi? There are, of course, many more questions along the same lines, 
e.g., what are probabilities of other more complex knots? what is the entropic response 
of a topologically constrained loop to various perturbations, etc. 

Most of what we know about these "probabilistic topology" questions is learned 
from computer simulations. In particular, it has been observed by many authors over 
the last 3 decades ^J El El H^l E! that the trivial knot probability depends on the 
length of the loop, decaying exponentially with the number of segments in the loop, 
N: 

Wt,i^ = Aexp{-N/No) . (1) 

For some lattice models this exponential law, in the iV — > oo asymptotics, was also 
mathematically proven [201 IS]- It was also noticed 17^ that the same exponential 
law, with the same decay parameter Nq, also describes the large N asymptotical tail of 
the abundance of any other particular knot - although for complex knots exponential 
decay starts only at sufficiently large N (as soon as the given knot can be identified 
as an underknot |24p. 

An alternative view of formula (^, useful in the context of thermodynamics, 
implies that the removal of all knots from the loop is associated with 
thermodynamically additive (linear in N) entropy loss of 1/Nq per segment; in other 
words, at the temperature T, untying all knots would require mechanical work of at 
least fcsT/A^o per segment. 

Another manifestation of the importance of the A'o parameter was found in 
the recent series of works |22| ESI ^M 1211- These works belong to the direction 
|26II27I 1^1291 I3()j addressing the spatial statistics of polymer loops restricted to remain 
in a certain topological knot state. It turns out that even for loops with no excluded 
volume and thus are not self-avoiding, A^o marks the crossover scale between mostly 
Gaussian {N < Nq) and significantly non-Gaussian {N > Nq) statistics. Indeed, at 
N < Nq, locking the loop in the state of an unknot excludes only a small domain 
of the conformational space which produces only marginal (albeit non-trivial [2H|) 



corrections to Gaussian statistics - for instance, mean-squared gyration radius of the 
loop is nearly linear in N. By contrast, at A^ > A^Oi the topological constraints are 
of paramount importance, making the loop statistics very much non-Gaussian, and 
consistent with effective self-avoidance 23, 24, 26, 30 . 

Thus, it seems likely that the parameter iVo might hold the key to the entire 
problem of knot entropy. We therefore decided to look at this parameter more closely 
in this paper. 

Present understanding of the values of Nq is quite modest. First, the constant's 
value was invariably found to be quite large, around 300 for all examined models 
of "thin" loops with no excluded volume, or no self-avoidance ^| El El 1221 1^ ■ 
Second, it is known that knots are dramatically suppressed for "thick" self-avoiding 
polymers, which means that TVq rapidly increases with the radius of self-avoidance 
[TtI IT!?] . The latter issue is also closely connected to the probabilities of knots in 
lattice models, where the non-zero effective self-avoidance parameter is automatically 
set by the lattice geometry. In the present paper, we will only consider the arguably 
more fundamental case of loops with no self-avoidance. 

The starting point of our analysis is the observation that Nq appears to be 
noticeably different for two standard polymer models for which common sense 
suggests that they should be equivalent. Both models can be called freely-jointed 
in the sense that they consist of N rigid segments with free rotation in the joints. 
However, in one model all segment vectors are of the same length, while in the other 
model segment vectors are taken from a Gaussian distribution. The motivation to 
consider the Gaussian distributed step vectors comes from the idea of decimation, 
or renormalization: we can start from the loop with Ng segments of equal length 
and then group them into A^ 3> 1 blobs oi g ^ 1 bare segments, each blob having 
nearly Gaussian distributed end-to-end vector. With respect to the knot abundance, 
the fixed length model was examined in El ^-nd the Gaussian model in El- It was 
noticed that Nq for the Gaussian distributed steps was larger than for identical steps, 
assuming no self-exclusion in both cases. No attention was paid to this observation, 
possibly because there was no confidence that the observed difference is real, in the 
context of the numerical error bars in the pertinent measurements. 

Recently, El EH) more detailed data became available which suggest that indeed 
A'^o is different for the two models, with fixed or Gaussian distributed steplength. A 
similar result was independently obtained by Vologodskii O]. Latter in this article, we 
present even better quality data supporting the same observation that Nq is different 
for these two models. This is a rather disturbing observation. Indeed, the idea of 
universality in polymer physics (S^ suggests that there should not be any difference 
between these two models as far as any macroscopic quantity is concerned. For 
instance, not only is the mean squared gyration radius the same for both models, 
but even the distribution of the gyration radii are the same, except far into the tails. 
In general, the difference between polymer models of this type becomes significant 
only in the strong stretching regime El or at high density El- Even if one takes into 
account the idea that knots, when present, are most likely localized along the chain 
El 1371 1381 1391 BH] , it is unclear how this fact can manifest itself for the loop that has 
no knots. 

Thinking generally about the loop models with fixed or Gaussian steplength, our 
reaction to this discrepancy is to realize that the major difference between the two 
freely-jointed loop models is that the Gaussian model may have a few unusually long 
segments, suppressing the ability of other shorter segments to wind around, and thus 
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Figure 1. (Color online) The fraction of loops of A'^ segments with trivial topology 
(Oi, or unknots) follows the decaying exponential form given in eq.0, as seen in 
this semi-logarithmic plot. The value of the decay constant A^o varies considerably 
as the loop structure is changed. This figure displays the fraction trivially knotted 
for loops of fixed or gaussian-distributed steplengths. The incongruity of A^o 
between these two models, which for a variety of metrics are indistinguishable 
at large loop lengths, is apparent. The figure also shows data for a modified 
random-flight model, eq. Q. In this model, the value of A^o is tunable by way 
of the parameter a, which allows a substantial range of A^o values, as indicated 
in the plot. The other model studied has bimodal distribution of steplengths, as 
described in section lij.3l . and although not shown in the figure, allows similar 
variability in the decay length A^o- 



decreasing the possibility for knots to occur, and accordingly, increasing iVo. Thus, the 
ability to take long strides might account for the comparative slowness of the ensemble 
of loops with Gaussian steplengths to diversify their knot spectrum with increasing 
A'^. The main goal of the present work is to investigate this conjecture. 

The plan of the work is as follows. After a brief description of our computational 
algorithms used to generate closed loops and to identify their topologies (section^ , we 
present computational results (section ^ on knot abundance for a variety of models 
differing in the width of their steplength distribution. In addition to the already 
mentioned loops with fixed, and Gaussian-distributed steplengths, in order to look at 
the even broader distributions which allow for very long segments, we also generated 
loops with the generalized Cauchy-Lorentz "random-flight" distribution. Finally, we 
include loops of bimodally distributed fixed steplength. 

In brief, our results are as follows. First, we confirm the exponential decay law 
of the unknot probability, formula (^, across all models examined. Second, we find 
qualitatively that indeed a wider distribution of the segment lengths leads to knot 
suppression, ie a larger TVg. Third, and most unexpectedly, we find that Nq does not 
show any signs of any singular behavior associated with the divergence of mean squared 
segment length or any other moment of the segment length distribution. Instead, Nq 
blows up and appears to grow without a bound when the distribution of segment 
lengths approaches the border of normalizability 



2. Models and Simulation Methods 

2.1. Models 

All polymer models referenced and employed in this work use a freely jointed model 
to represent a polymer loop. The polymer is represented by a set of N vertices in 
3D, with position vectors Xi, where the step between successive vertices is described, 
f = Xi+i — Xi. In all models, we assume that the distribution of segment vectors 
r, which we call P{r), is spherically symmetric and depends only on the steplength 
r = |f|, such that 

{r) ^ fP{r)<fr = . (2) 

We also assume that the mean squared steplength is always the same (when defined! 
- see below), we denote it i: 

(r^) = I r^P{r)cPr ^ in f r*P{r)dr = i^ . (3) 

With this in mind, the simplest measure of the distribution breadth involves higher 
order moments: 

^J^'>-'f> . (4) 

where (r^) = 47r/„°° r^P[r)dr. 

Specifically, we analyzed the following models. 

The fixed steplength model is described by the distribution 

For this model, of course a — Q. With N segments, the loop's contour length is 
obviously L = Ni and the mean squared gyration radius of the loop is {R'i) — 
{N + 1)^/12. 

The gaussian steplength model is generated by the distribution 
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in this case, a = ■\/2/3. The contour length of the A^-segment loop in this model 
is L — N {If]) — N£y^8/37r and the mean squared gyration radius is (R'i) = 
{N -1/N)e/12. 

The random flight steplength model is obtained from the generalized Cauchy- 
Lorenz distribution (also known in the theory of Levy flights 4S^) of the form 

asinf37r/a]/(4^3^2\ 
Pir] = — ^ '- f7) 

where the factor in the numerator ensures normalization. Here, a, a parameter of 
the distribution, must be greater than a > 3, otherwise the normalization integral 
diverges. Nevertheless, it would be fair to speak about a family of random-flight 
models, parameterized by a > 3, instead of just one model. Varying a allows us to 
work with a "tunable" distribution. These distributions' "fat" power-law tails lead to 
diverging moments (which is why they are used to describe super-diffusive behavior 
seen in biological foraging 42i i42i iMi dSl EH EH j and quite recently, in the diffusion of 



bank-notes across the United States, @B1)- Specifically, the ensemble-averaged contour 
length of the loop is well defined only at a > 4 {L — Nlsin {Sir /a)/ sin (47ra)), mean 
squared gyration radius exists at a > 5, and a only exists at a > 7, in which case it 
is equal to 



sin f57r/al , , „, ,„, 

^ ' ^ 1 (at a > 7) . (8) 



y sin [Stt/q:] sin [77r/Q!] 

Finally, we also include loops with bimodally distributed steplength. For 
these loops, 

which means that two possible steplengths, £i or £2, occur with probabilities Pi and P2, 
subject to the normalization conditions, P1 + P2 = 1, and (r^) = Pi^f -1-^2^2 = ^^- ^11 
bimodally distributed models can be conveniently parameterized by Pi and A = ^2/^1- 
For these models. 



P^ + ('-^^^'' -1, (10) 



y (Pi + (1 - Pi)A2)2 

might be very large if A is very large and Pi is rather close to unity (A^ » 1/(1 — Pi) ^ 
!)• 

2.2. Loop generation 

Unbiased generation of closed loops is of decisive importance for our work. Recently, 
we gave a detailed review of the existing computational methods to generate 
statistically representative closed loops (see the last section of the work [23 )■ In 
principle, the best way to generate loops is based on the so-called conditional 
probability method. The idea is that a closed path is generated as a random walk, 
step by step, except after the completion of k steps, the next step, fc -I- 1, is generated 
from the analytically computed probability distribution of the step vector r, subject 
to the condition that after N — k more steps, the walker returns to the starting point. 
This idea was first suggested and implemented for Gaussian distributed steps |1T] . 
Recently, we implemented this method for steps of equal length j25| . Unfortunately, 
this method is computationally costly, and appears to be prohibitively difficult to 
implement for more sophisticated models, such as random flight. 

We therefore use the simpler method, called the method of triangles. This method 
123 ^M Uni generates loops of N segments with A^ divisible by 3. It involves creating 
a set of 7V/3 equilateral triangles, each randomly oriented in 3D space. Each triangle 
is considered a triplet of vectors with zero sum. A random permutation of the N 
edge vectors which make up these A^/3 triangles, and then connecting all N vectors 
head-to-tail, creates a loop which will be closed, as the N bond vectors together have 
vector sum. Of course, this method imposes correlations between segments. We 
therefore take special care to compare the results of this method with the unbiased 
generation using the conditional probability method for both Gaussian distributed and 
equal length step models. We found that no appreciable deviations in knot abundance 
data arise from the imperfection of the triangle method. We therefore use the method 
of triangles to generate the random flight and bimodal distributed loops, for which no 
alternative method is available. To avoid even the slightest problems with correlations. 



implicit in our triplet method, and to ensure that the decay of trivial knot probability 
is in the exponential regime, we exclude the data from small loops and fit the trivial 
knot probability on the interval N £ [50, 300]. 

2. 3. Identification of topology and statistics 

For each of the models listed above, we generated loops of up to at least N = 300 
segments. 

Once the loop was generated, its knotted state was assessed computationally with 
the Alexander determinant, | A(— 1)|, as well as the Vassilev invariants of degree 2 and 
3, v2 and u3. For this purpose, we employed knot analyzing routines described in 
detail elsewhere I33| . 

The fraction of generated loops with trivial topology was recorded for each loop 
type. In the interval 50 < A^ < 300, every sample consisted of at least 10^ loops. As a 
result, the plot of trivial knot probability was created for each model, with statistical 
error bars smaller than the data points in Figure ^ Based on the data, A^o was 
measured for every model in our repertoire. 

3. Results 

Our main results are shown in figure ^ There, we present the semi-log plots of the 
data on the trivial knot probability as a function of the number of segments N for a 
variety of models. To begin with, all our data agree quite well with the exponential 
character of the wtrwiN) dependence, formula (^. Our main emphasis is therefore 
the study of the characteristic value of A^o for various models. 

3.1. Loops of fixed steplength versus Gaussian distributed steplength 

Let us start with the two most commonly used models. In the case of fixed length 
steps, we obtain A^o ~ 240 (see also US]). By contrast, for Gaussian distributed length 
steps we get A^o ~ 305 (see [2H1)- 

All data for the Gaussian model was generated using the fundamentally unbiased 
conditional probability method. For the model with fixed length steps, we compare 
in figure |21 the data obtained by the conditional probability method "^ and data 
generated by the much more efficient method of triangles (see above section IT^ . As 
the figure indicates, there is practically no visible difference in the results. Accordingly, 
we unreservedly rely on the triangles method in the rest of this work. 

Although our main attention in this paper is on the trivial knots (unknots), we 
show in figure 13 some data for more complex knot probabilities. Our data at least 
do not contradict the assertion that the probability of every particular knot decays 
exponentially at sufficiently large A^, with the same decay length as the trivial knot 
probability. At the same time, our data also confirm the systematic difference between 
models - decay length, although the same for all types of knots, does depend on the 
distribution of step lengths involved. Specifically, probabilities of knots 3i and 4i 
decay at a noticeably slower rate for the Gaussian model than for the model with 
fixed steps. 
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Figure 2. (Color Online) Semi-log plots of trivial knot probability against the 
number of segments for the Gaussian distributed step length model (□), and for 
the fixed steplength model, the latter generated either via conditional probability 
method (o) or triangles method (A). The results indicate no dependence on the 
generation method, thus lending credence to the simpler method of triangles. The 
results also indicate an almost perfect fit to the exponential law, Eq. 0, with 
A'^0 ^ 240 for fixed steplength and TVq ~ 305 for Gaussian distributed steps. 




Figure 3. The probabilities of simple non-trivial knots 3i and 4i are plotted 
against the number of segments in the loop in semi-log scale. The data at large TV 
is consistent with the exponential form of these probabilities (plotted as solid lines 
in the figure) , which are characterized by the same parameter A^o found for trivial 
knots. As with the case of trivial knots, we see the difference between the fixed 
steplength model and Gaussian distributed model, namely, that any particular 
knot probability decays slower for Gaussian distributed segments than for fixed 
segments. 



3. 2. Random flight loops 

The results of our study of loops with random-flight steplength, a E [3.5,30], are 
summarized in figures ^ and ^ The raw curves of probability given in figure ^ 
clearly show that the odds of finding an unknot in a set of loops get increasingly 
unfavorable as a decreases. That is not unexpected: at smaller a, the probability 
distribution JTJ) acquires an increasingly fat tail, which implies the presence of a 
fraction of exceptionally long segments, and they of course suppress the chance of 
knots. 

At very large a, the knot probability for the random flight model appears similar 
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Figure 4. (Color Online) Values of the trivial knot decay length, A^o, measured 
from loops constructed of steps with generalized Cauchy-Lorenz distribution 
(reminiscent of Levy-flights). The decay constant is expressed as dependent on 
the variable a which defines how fat the tail of the distribution remains at large 
steplengths. The knot probability decay constant, A^o, approaches a constant 
value of about 266 at large a. The inset shows the same data in double logarithmic 
scale. It is seen that A^o blows up approximately as (a — 3)~^ when a decreases. 
It is interesting to note that the dependence of A^o on a does not show any signs 
of irregularity as a crosses values at which various moments of segments length 
distribution start diverging (for instance, mean squared gyration radius diverges 
at o < 5 and mean contour length diverges at a < 4). 



to the data for fixed steplength loops. Indeed, as figure ^indicates, A^o for the random 
flight model at very large a approaches iVo ~ 266 which is not dramatically different 
from A^o ~ 240 for the fixed length steps. In fact, the remaining difference might be 
associated with the fact that even at very large a, the random flight model, although 
it has essentially no very long steps, has some relatively short ones, which might 
account for the discrepancy in iVp. Figures ^ and 01 show further that at a « 7.5, the 
random- flight loops behave essentially the same way as Gaussian steps in terms of TVq. 
Most interestingly, figure 01 shows no sign of anything unusual happening to A*"o 
at the values of a at which various physically important moments of the segment 
length distribution {T)) start diverging. For instance, at a < 5 the mean squared 
gyration radius diverges, at a < 4 even the contour length of the loop diverges - and 
yet none of these facts find any visible reflection on the dependence of A^o on a. A^o 
keeps smoothly increasing with a, with a maximum measured value of A^o ~ 4800 at 
a = 3.5. It appears that A^o in fact blows up and goes to infinity as a approaches 
3 - the border below which the distribution is not normalizeable. Moreover, as 
the inset of figure ^ shows, this divergence is well approximated by the power law 
dependence of the form, 

N^^Nl^ Bia-a*)-^ , where 

A^o « 266 , B « 1080 ± 150 , a* = 3 , and 

/3 « 2.1 ±0.1. (11) 

In fitting the data with this power law, we ignored the small irregularities visible 
around a « 13 (or In [a — ao] ~ 2.2), which we attribute to the numerical problems 
with the cutoff implemented in our use of eq. (jA.l|l . 

The meaning and physical origin of the apparent criticality observed at a 
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Figure 5. (Color Online) The surface plot of A^o measured for loops of bimodally 
distributed steplength. The two steplengths used, £i and ^2, related by A = £2/^1, 
occur with probabilities Pi and P2. The surface is described by two degrees of 
freedom, A and Pi . The model reduces to the fixed steplength model of section 
I3.1l in the conditions Pi = 0, or Pi = 1, or A = 1, as is seen by the reduction of 
A^O ~> 241 in these limits. Interestingly, the bimodal distribution of steplengths 
never reduces A^o below the limit of No = 241 of the fixed-steplength model. 
The maximum value of A^o which we were able to observe, A^o = 974, occurs at 
Pi = 0.15, A = 0.01. 



approaching its minimal possible value of 3 currently evades our understanding. 



3.3. Loops of bimodal- distributed steplength 

For this model, we examined the interval N S [50, 600], with at least 10^ loops in each 
simulation record for N < 300, and at least 10^ loops in those records used with length 
N > 300. We also consider parameters Pi and A = -^1/^2 in the intervals Pi S [0, 1] 
and A G [0.01, 1.0], respectively. By symmetry, it is sufficient to look at A < 1; this 
means. Pi is the fraction of shorter segments. 

The raw data from the loops with bimodal steplength are presented in figure El 
This surface plot charts the change in Nq as a function of two parameters of the model. 
Pi and A. Below the surface is the corresponding contour plot of Nq as a function of 
these two parameters. The changes in Nq are smooth, and there do not seem to be 
any singularities in the behavior of A*o. The maximum observed value of Nq = 974 
occurs at Pi = 0.15, A = 0.01. This maximum of Nq appears to be rather sharp, small 
deviations in Pi from 0.15 lead to smaller Nq values. From the data, we believe that in 
this model A'o is maximized when the fraction of shorter segments is Pi — 0.15 ± 0.05. 
As regards the second degree of freedom. A, while our data certainly shows that as 
A — * 1 the bimodal system approaches the fixed steplength system of section 13.11 (ie 
A^o = 241), in the opposite direction of A ^ 0, it is not clear if Nq will continue to 
increase in an unbounded way or be in some way encumbered. In this regard our 
present knot analysis machinery is limited by the relative disparity between segments 
of different length, and more work needs to be done to elucidate the scaling of Nq in 
the limit A ^ 0. 

Qualitatively, all of our data are consistent with the idea that what suppresses 
knots is the presence of a fraction of unusually long segments. One could then 
hypothesize that Nq might depend on some unique property of the segment length 
distribution, for instance, a, as defined in eq. Q. This hypothesis is tested in figure 
El the figure indicates that the hypothesis fails. Nevertheless, the results presented 




Figure 6. (Color Online) The value of A^q for several different loop models is 
displayed against the steplength distribution width a, which is defined in eq. l4l . 
The data in this figure fail the hypothesis that a larger a universally maps to a 
larger A^o- The measurement for fixed steplengths, shown in the lower left corner 
is the point from which all bimodal loop data emanates. As labeled in the figure, 
the manifold for bimodal steplength is indexed by lines of constant Pi and A. 



in this figure are interesting, as they show that large A^o can be achieved by the 
combination of segment length difference (A) and segment type fractions. 

Although a functional relationship between A^o and a seems evident in our data for 
random flight loops, the data we have for loops with bimodal distribution of steplength 
do not suggest a simple single parameter which determines A^o- 

4. Conclusions 



In a qualitative sense, it does seem that there exists a relationship between A^o ^^nd 
the reach of successive segments within the chain (as seen in Figure El ■ It seems 
qualitatively clear indeed that knottedness is greatly suppressed by the presence 
of some very long segments. Thus, the slowly decaying tail of the segment length 
probability distribution, or the presence of a small fraction of very long segments, 
implies a large Nq, which sounds natural. However, our understanding of these 
observations beyond the qualitative level is limited. We found that A^o appears to 
exhibit no singularity associated with divergence of any natural characteristics of the 
loop, such as its gyration radius or even contour length; instead, A^o exhibits power 
law critical behavior when the segment length distribution approaches the boundary 
of normalizability. In addition, we do not know which property of the segment 
length distribution determines A^o- We have established that this is not simply the 
distribution width, and from the results with random flight distribution it seems also 
clear that it is not based on any finite moments of the distribution. We consider the 
development of a fuller understanding of the variance of Aq a compelling challenge. 

The data clearly show that wide variation in the behavior of A^o is possible, 
including very large values of A^o in certain models. Given that A^o plays the role of 
the cross-over length for the critical behavior of topologically constrained loops |2SI j 
we can speculate that the models with very large A^o are in some way similar to the 
models of self- avoidance in the vicinity of the ^-point. We think that this analogy 



deserves very close attention. 
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Appendix A. On the numerical implementation of generalized 
Cauchy-Lorenz distribution 

To generate steplengths from the random-flight distribution , we first take a random 
number q from the uniform distribution on the interval [0, 1] and then find steplength 
r as r = f{q), where the mapping f{q) is determined by the equation 

rfil) 
q= P{r)ATrr^dr , (A.l) 

Jo 

where P{r) is given by eq. Q). Although a closed- form representation of the right- 
hand side of eq. (|A.1|I exists in the form of an Incomplete Beta Function, [201, we 
chose to implement the mapping via a numerical interpolation of tabulated values of 
the integral. As the power-law tail of this distribution is quite fat, particularly at small 
a, accurate representation of the integral becomes challenging at large steplengths. 
This work is ultimately numerical, and we are forced to truncate the representation 
of the integral to 1.0 at some maximum steplength, ie specifying an upper bound on 
f{q)- The slow convergence of the tail is behind the appearance of the few outliers in 
the data for the random- flight method, figure^ 
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